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ABSTRACT 

We study the stellar response in a spectrum of potentials describing the barred spiral galaxy 
NGC 1300. These potentials have been presented in a previous paper and correspond to three 
different assumptions as regards the geometry of the galaxy. For each potential we consider 
a wide range of Q^, pattern speed values. Our goal is to discover the geometries and the Q^, 
supporting specific morphological features of NGC 1300. For this purpose we use the method 
of response models. In order to compare the images of NGC 1300 with the density maps of 
our models, we define a new index which is a generalization of the Hausdorff distance. This 
index helps us to find out quantitatively which cases reproduce specific features of NGC 1300 
in an objective way. Furthermore, we construct alternative models following a Schwarzschild 
type technique. By this method we vary the weights of the various energy levels, and thus the 
orbital contribution of each energy, in order to minimize the differences between the response 
density and that deduced from the surface density of the galaxy, under certain assumptions. 
We find that the models corresponding to Q.p ^ 16km s ' kpc ' and Qp w 22km s ' kpc ' 
are able to reproduce efficiently certain morphological features of NGC 1300, with each one 
having its advantages and drawbacks. 

Key words: Galaxies: kinematics and dynamics - Galaxies: spiral - Galaxies: structure - 
ISM:kinematics and dynamics 
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1 INTRODUCTION 

In iKalapotharakos et al.l jlOlOl) (hereafter Paper I), we have pre- 
sented three different general models representing the potential 
of NGC 1300. The morphology we investigate can be observed 
at a deprojected near-infrared image of the galaxy (Fig. [T). Each 
of these model-potentials corresponds to different assumptions re- 
garding the distribution of the luminous matter of the galaxy in the 
third dimension, perpendicular to the equatorial plane. They repre- 
sent limiting cases for the geometry of the system, and vary from 
the pure 2D (Model A) to 3D cases. The two 3D cases correspond 
either to a pure cylindrical geometry of a 3D disc with a constant 
scale height (Model B) or to a combination of a spheroidal com- 
ponent representing the central part of the bar with the cylindrical 
geometry of a 3D disc for the rest of the luminous mass (Model C). 
All of them have the option of the inclusion of two additional terms 
in the potential representing the central mass concentration and the 
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dark halo component respectiv ely. These term s are c onstrained by 
the kinematical data derived bv lLindbladetal.l ( ll997l) . 

We note that in our calculation we do not take into account 
the outer spiral arms that exist beyond the edges of the frame of 
Fig-in These spiral extensions are very weak in the near-infrared, 
while they are conspicuous in the optical. This indicates that they 
consist mainly of young objects, thus their contribution to the mass 
distribution of the galaxy is small (Grosb0l et al. in preparation). 

In the present paper we investigate the detailed dynamics in 
the three cases using response models (hereafter RM). As we have 
already mentioned in Paper I the global dynamical behavior in ro- 
tating galaxies, is crucially determined by the assumed value of 
the pattern speed SI,,. Thus, in this paper we use the method of 
response models (Patsis 2006) in order to find which CI,, values, in 
each Model (A, B and C), are able to reproduce the various morpho- 
logical features of NGC 1300. We construct our models under the 
assumption of a single pattern speed and we will assess the results 
at the end. For this reason, we introduce an index that is a gener- 
alization of the Hausdorff distance (see e.g. Deza & Deza 2009) so 
that we can quantify the (dis)similarity between the density maps 
of the response models and the K-band deprojected NGC 1300 im- 
age. Comparisons that involve the surface density of the galactic 
disc are done under the assumption of a constant M/L ratio (see 
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Figure 1. The neai'- infrared surface brightness of NGC 1 300, in gray scale 
(color scale in online version). This image has been obtained after applying 
all the necessary amendments (see Paper I) and after deprojecting the image 
with the considered position and inclination angles (PA, IA)=(87°, 35°). 
The solid curves ai'e six iso-density contours describing the basic contour 
shapes of NGC 1300. The adopted distance is D=19.6 kpc. 



Paper I). We remind also, that from the image we used (Fig.[Tll we 
have removed point-like sources, that are likely to be young stellar 
clusters (Paper I). 

This study allows also the comparison between various re- 
sponse models and consequently indicates, which kind of geome- 
try applies better to the NGC 1300 case. Nevertheless, our ultimate 
goal is the detection of the dynamical mechanisms behind the ob- 
served structures (bar and spiral arms). This needs a detailed anal- 
ysis of the orbital stellar dynamics for each case and is presented in 
IPatsis et all 1 I2OIOI) (hereafter Paper III). 

Our paper is structured as follows: in Section[2]we present the 
method of response models and the corresponding results. In Sec- 
tion|3]we describe the generalized Hausdorff distance that helps us 
quantify the similarity between the various response models and 
the galaxy. In Section|4]we check the ability of the response mod- 
els to describe better the morphological features of NGC 1300 by 
varying the relative contribution of the energy levels. Finally, we 
discuss our conclusions in Section[5] 



2 RESPONSE MODELS 

These models show the response of an initially axisymmetric stellar 
disc, with particles moving at the beginning of the simulation in 
circular orbits determined in the axisymmetric part of the potential. 
The steps we follow in our numerical experiments are: 

(i) we choose a potential (among Models A, B, C) and the pat- 
tern speed value dp, 

(ii) we populate uniformly the disc up to i?„,„, = 15kpc (the 
(x,,y,) of 10* particle positions on the plane are taken at random 
positions). 
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Figure 2. The RcR/Rtar ratio as a function of the pattern speed for our 
models A (filled symbols) and B (empty symbols). Squares (blue colored in 
the online version) and circles (red colored in the online version) correspond 
to the long and the short semi-major axis of the bar, respectively. Note that 
there are ranges of Clp values that give multiple Lagrangian points. 



(iii) we set each particle in circular motion (in the rotating 
frame) with velocity Vdn- 
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Clp R, where Oq is the 

'K K 

axisymmetric part of the potential, 

(iv) we grow the non-axisymmetric terms of the potential lin- 
early from towards their full amplitude within two pattern rota- 
tions, 

(v) we integrate the particles' orbits for 25 pattern rotations (we 
consider that a particle has escaped and we stop the integration of 
its orbit when it reaches a radius R > 22 kpc). 

(vi) Finally we construct density maps, by converting our data 
files to images. For this we consider a grid and we take into account 
the numerical density of the test particles on the disc. 

In Paper I we have studied the distribution of the Lagrangian 
points in all general models (A, B, C) and for an extended range 
of Q.,, values (see fig. 10 of Paper I). In Fig. [2] we plot the ratio 
RcRlRbar of the corotation radius Rcr over the semi-major axis of 
the bar Ri,„r as a function of the corresponding Q.p value. We have 
omitted the points corresponding to Model C since they are very 
close to those of Model B. We have considered the bar extending 
up to the outer white isocontour of Fig.[T](third white contour start- 
ing counting from inside). This contour is not symmetric to the 
origin, since the right semi-major axis of the bar is longer than the 
left one. Thus, we set Ri,ar = 8.9 kpc for the right semi-major axis 
of the bar and = 8.0 kpc for the left one. From the above it 
is obvious that the ratio RcRlRbar depends on the side of the bar. 
The considered Rqr value is the radius of the Lagrangian points 
(Li,L2) lying in the corresponding side. In Fig.[2]the filled circles 
and squares correspond to Model A, while the empty ones to Model 
B. The squares (colored blue in the online version) and the circles 
(colored red in the online version) represent the ratio RcRlRbar on 
the right and on the left side (Li and L2 area), respectively. It is 
argued that Rtar should not be longer than the corotation radius 
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Figure 3. The surface density in gray scale (color scale in the online version) of the RMs coiresponding to the potential of the pure 2D case (Model A). At 
the top of each panel is indicated the Clp value. In each frame we plot over the density maps of the RMs six representative contours (solid lines) of the K-band 
image of NGC 1300 (the same with those plotted in Fig.[T). The density scale is the same as in Fig.[T] 
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Figure 4. As Fig. [3] but for the potential con'esponding to the cyhndrical geometry (Model B). 
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RcR- This argument is sustained by dynamical studies according to 
which the orb ital content beyond corotation is not able to support 
the bar shape t Contopoulod l 1980). In our models, as we have men- 
tioned already in Paper I, we see that there are ranges of £1,, values 
with multiple Lagrangian points (especially in Model A). Particu- 
larly, in some cases the ratios RcRlRtm corresponding to multiple 
Lagrangian points are quite different when we consider as corota- 
tion distance the distance of the one or the other Lagrangian point. 
E.g. in Model A at fl « 25 km s"' kpc"' the ratio RcR/Rbar varies 
from X 0.5 (unphysical value since it is significantly lower than 1) 
up to » 1.1 (reasonable value). These cases are quite complicated 
and lie beyond the standard paradigm. Thus, there is no a priori 
knowledge about the orbital behavior and the morphologies they 
can support. In any case, the models with Clp > 28 km s"' kpc"' 
have single Lagrangian points corresponding to significantly low 
RcRlRbar values. 

In this study we present the results corresponding to the Q,, 
values from 12 to 26 km s"' kpc"' . Empirically it has been also 
realized that beyond the limits of this range the response mod- 
els were particularly problematic in the sense that there was ob- 
viously no good agreement with the image of the galaxy. The mod- 
els with Sip > 17km s"' kpc"' correspon d to "fast bars" as are 
usually called these with Rc RlRtm < 1.4 jRautiainen et al.|[2008l : 
iDebattista & Sellwoodl200(j) . 

In Fig.[3]we plot, in gray scale (color scale in online version), 
the density maps we obtain from the RMs corresponding to the 
potential of Model A (2D case) for the 12 < Hp < 26 km s"' kpc"' 
values. Here, we plot the RMs corresponding to the 8 indicated Q,, 
values. However, we have data for models with totally 15 Q,, values 
within the above range of Q.p. On each panel of this figure we have 
also overplotted the same (as in Fig. [TJ six iso-density contours 
(solid lines) of the K-band image of NGC 1300, outlining its major 
morphological features. 

In this figure we observe that the low fi,, values (first row) 
fail to reproduce both the bar and the spiral structure. The response 
bar is conspicuously broader than the galactic bar, while the spi- 
ral arms are absent or flocculent. For higher values of Sip we get 
spiral arms in good agreement with the real ones (see panel for 
Sip = 16 km s ' kpc"' ) although the bar is still broad. For even 
higher Sip values (third row) the formed spiral arms of the models 
do not match with the spirals of the galaxy. Nevertheless the bar 
clearly gets an ansae character in agreement with the NGC 1300 
bar morphology (see panel for Sip = 22 km s"' kpc"' as best ex- 
ample). Finally, the highest Sip values (last row) produce very open 
spiral arms and/or rings and a disintegrated bar. 

Figure|4]presents the RMs for Model B (thick disc, cylindrical 
geometry) at the same Sip values as for Model A (Fig. [3}. In this 
case the response bars have always high density regions at their 
ends, underlying a clear ansae-type character. A common feature 
is also a density enhancement along the minor axis. The bar mor- 
phology has a good relation with the galactic bar for low Sip val- 
ues (two first rows). This indicates that a proper selection of the 
bar- supporting orbits might lead to a best fitting of the NGC 1300 
bar. As Sip increases (panels in the third and forth rows) the re- 
sponse bar becomes shorter than that of the galaxy. The response 
spiral arms are not always well described but there is a range of 
Sip = 16-22 km s"' kpc"' , where occasionally we find high 
densities in the models at the regions corresponding to the galactic 
spiral arms. 

Figure|5]is similar to Figs.[3]and|4]but now for Model C (spher- 
ical and cylindrical geometry). The morphology of the response bar 
resembles more the morphology we get in Model A especially for 



the lower SI,, values (first two rows). However, the inner structure 
of the NGC 1300 bar {R < 2kpc) is fitted clearly better in this case. 
The general spiral response is closer to that of Model B while it is 
in general fainter. 

Up to this point, the similarity between the various RMs and 
the morphological features of the galaxy have been made by eye. 
Although the eye is usually a good selecting tool it lacks of objec- 
tivity and capability of quantification. For this reason, in the next 
section, we introduce an index in order to quantify the comparison 
between the NGC 1300 morphology and the response models. 



3 QUANTITATIVE COMPARISON 

In Fig. [T] we have plotted six iso-density contours of NGC 1300 
reproducing its basic morphological features. Our goal is to "mea- 
sure" in an objective way the similarity of an RM density map with 
the galactic morphology. Thus, we want to establish a quantitative 
criterion for this resemblance. 

The methodology of this study is related with image compar- 
ison and pattern recognition techniques. A simple and well known 
index used in such studies is the "Hausdorff Distance" d/^, that is 
defined as follows: let A, B be two non-empty subsets of a metric 
space (M, d). In this case the Hausdorff distance d,i between A and 
B reads 



dn = max sup inf d{a, b), sup inf d{a, b) 



(1) 



Considering two sets of points corresponding to two contour lines 
and d being the Euclidean distance, measures the maximum of 
the distances between the points of each curve from the other one 
(as a whole). It is obvious from the above definition that dn = 
corresponds to identical lines. The major issue about rf^ is its sen- 
sitivity to outliers. Let us suppose that a contour line consists of two 
parts, the major one being identical to a contour line of NGC 1300, 
while the minor one being a tiny part far away. In such a case al- 
though we have a good matching of the response model with the 
galaxy morphology (at least for the specific contour), dn reads a 
false high value. There have been proposed many modifications of 
du that aim to reduce this drawback (Zhao et al. 2005, and refer- 
ences therein). Here we introduce an index which is a generaliza- 
tion of d^ and it is appropriate for our study. 

We call this index "Generalized Hausdorff Distance", dgu, and 
it is defined as follows: let Ci, C2 be two curves on the plane of non- 
zero length and > respectively. The dcH between Ci and C2 
is 



inf d(a, b) 



-dl + 



inf d{a, b) 



dt, (2) 



where the integrations are over the curves Ci and C2 and |HI denotes 
the Euclidean norm. The integration adds the contributions from all 
points, while the division with the curve lengths and the norm 
aims to the normalization of the final doH value so that it can be 
comparable for two totally different contour lines. In practice, for 
the computation of doH we sample both curves with a considerable 
number of points Nc^ , Nc-, and we calculate the quantity 

Nc. 



NC2 

What we do is to find which contour line of each RM corresponds 



1 1 1 1 

doH = / — mm d(ai,b)-\ > — mind(a,, i). (3) 

beCi 
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Figure 6. In the left-hand column we plot for the RM, the minimum generalized Hausdorff distance dan value relative to the NGC 1300 contour line shown in 
the right-hand column (same row). In the middle column we plot the minimum dan values corresponding to SMs (see Sect.|4). Note that filled circles, empty 
circles and empty squares (red, blue and green colored lines in the online version) correspond to pure 2D geometry (Model A), to cylindrical geometry (Model 
B) and to spheroidal+cylindrical geometry (Model C), respectively. The star symbol (light blue line in the online version) in the left-hand column coiresponds 
to a scenario with time-dependent Hp (see Sect.|5). 
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Figure 7. The six representative contour lines of NGC 1 300 (solid lines) together with the six contour lines con'esponding to the minimum doH values (dashed 
lines) for all the SMs of Model A. The comparable lines have the same color There are cases that these models are better than RMs in reproducing some 
morphological features of NGC 1 300. 



to the minimum dan value relative to each one of the contour lines 
of NGC 1300 we have plotted in Figs. [TT3l5] 

In Fig.[6](left-hand column) we have plotted the minimum 
values corresponding to the NGC 1300 contour lines shown in the 
right-hand column of the same figure. Filled circles, empty circles 
and empty squares (red, blue and green colored lines in the online 
version) correspond to Models A, B and C, respectively. Note that 
in this figure we have also included the models corresponding to 
Sip values not plotted in Figs. l3l5l Based on this criterion, the basic 
trends that are concluded by observing the behavior of dcH for all 
the Models (A, B, C) and for the studied range of Q,, values are: 

• the lower Q.p values in Model A are more efficient (than the 
higher ones) in portraying the outer contours corresponding mainly 
to the spiral structure of NGC 1 300 (see the two first rows of Fig.[6]l 

• the higher Q,, values in Model A are more efficient (than the 
lower ones) in portraying the inner contours corresponding to the 
barred structure of NGC 1300 (see third to fifth rows of Fig.[6]( 

• Model B is better than Model A in portraying the spiral struc- 
ture in high Q.p values, while for low Qp's, Model A is in general 
better than B, especially at the outermost isocontour (first row of 
Fig.|6j. 



• Model C is in general worst than Models A and B in portray- 
ing the spiral and bar structure, while it is definitely superior in 
portraying the very inner area especially for the middle and high 
Q.p values (see last row of Fig.|6f. 

Nevertheless, we should have in mind, that the quantification of the 
similarity of the contours of the model with those of the galactic 
image by means of dan, cannot in general account for small mor- 
phological differences, essentially in the phases, of the two curves. 
Despite this drawback dati is a reliable index in most cases encoun- 
tered in the present study. 

All the above results are related with the response models and 
depend to a large extent on their initial set up. These models provide 
useful information about the underlying dynamics and the morpho- 
logical features they can support. Despite the fact that the setting 
up of the models is done in an unbiased way, there is no physical 
reason to exclude weighting some energies more than others in an 
effort to improve the similarity between RMs and galaxy morphol- 
ogy. In order to explore this we modify our models by introducing 
weights attached to each energy level. The energies of the particles 
are determined by giving them circular velocities in the axisymmet- 
ric potential. This determines where they will be end up in the well 
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Figure 8. Similar to Fig.|7]but for the potential corresponding to the cylindrical geometry (Model B). 



of the full potential. However, one can consider a totally different 
energy distribution of the particles of the model. The procedure and 
the corresponding results are described in the following section. 



4 WEIGHTING ENERGY LEVELS 



ISchwarzschilj ( 1 19791) in a pioneer study, presented a method 
for constructing time independent self-consistent solutions 
of sp ecific models. H i s met h od has been widely used since 
then jRichstond ll98Cl 1 19821: ISchwarzschildl 1 19821: iRichstond 
1984 iR ichstone & Tremaine' 'l 984l: ISch warzschild l Il993l: 
Merritt & Fridman 1996; Cretto netall 1 1999: Anto nini et al.l 



20091) mostly for models representing elliptical galaxies. The 



basic idea of this method is very simple: One considers a density 
distribution function and then finds the associated potential via the 
Poisson equation. The next step is to calculate a library of orbits 
corresponding to this potential for the desired energy range. By 
assigning a weight to each orbit one tries with the superposition 
of the weighted orbits to reproduce the initially imposed density. 
For models of nor mal spiral galaxies , a similar method has been 
used in the past by IPatsis et all ( Il99ll) for the investigation of the 



self-consistency in 12 galaxies, while in barred-sp iral systems has 
been applied bv lKaufmann & Contopoulo3 ( Il996l) . 

In practice what one is looking for is a way to minimize the 
difference between the imposed density and that corresponding to 
the superposed orbits. This detailed task is very demanding and is 
beyond the scope of this paper. Here we will construct our "best" 
models acting in a more rough way. We follow the rationale of 
Schwarzschild's method but we assign weights to the various en- 
ergy levels. The density distribution for each energy level is con- 
sidered being the one we get from the corresponding RM. This as- 
sumption is true only for phase space domains where Chaos dom- 
inates and for time scales for which the corresponding initial con- 
ditions have been evolved and are close to a dynamical equilibrium 
due to chaotic mixing. For phase space domains (of specific en- 
ergy levels) with significant fractions of regular orbits the above 
assumption is just a simplification since the distribution of the par- 
ticles on the existing tori can be supposed to be different than that 
corresponding to the RMs. 

The determination of the weights is made by the following 
procedure. Let Ng be the number of Jacobi constant (energy) inter- 
vals, which we give to the particles of a response model and the 
number of grid cells we divide the plane (x,y). In this situation we 
derive the weights vv,, (/ = 1, . . . , A',, w(i) > 0), that minimize the 
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Figure 9. Similai' to Fig.|7]but for the potential corresponding to the (spheroidal + cylindrical) geometry (Model C). 



quantity 

j=\ V ,=1 / 

where nig., is the mass contribution of the ;-th energy interval on 
the j-th grid cell and nij the mass contribution we consider corre- 
sponding to the same region at the NGC 1300 image on the j-\h 
grid cell. We consider the ratio me../mj in order to treat equally 
high and low density regions. The final considered surface density 
distribution is given by 

di=^ j=l,...,M. (5) 

where dj, S j are the density and the surface area of the j-th grid 
cell, respectively. In the present study, we used = 30 and A'^ = 
24 X 24 = 576 for the central 576 kpc^ area. 

Thereupon, we repeat the procedure of the determination of 
the contour lines of each of these "Schwarzschild type" models 
(hereafter SM) corresponding to the minimum cIgh value relative 
to each contour fine of NGC 1300 plotted in Figs.[T]and|3l5l 

Going back to Fig. [6] in the middle column, we plot the min- 
imum dan values, similar to the panels in the left-hand column 



(RMs), but for the density distribution corresponding to Eq. ^ 
(SMs). The comparison between the RMs and SMs behavior of 
dcH, shows the following: 

• In 2D geometry (Model A, filled circles or red color line in the 
online version) the SMs with the high Qp values are more efficient 
(than the corresponding RMs) in portraying the spiral structure of 
NGC 1300 (panels in the two first rows in Fig.[6j. 

• In cylindrical geometry (Model B, empty circles or blue color 
line in the online version) the SMs with the low Qp values are 
more efficient (than the corresponding RMs) in portraying the bar 
of NGC 1300 (third to sixth rows of Fig.|6). 

• In spheroidal+cylindrical geometry (Model C, empty squares 
or green color line in the online version) the SMs with the lowest 
Q.p values (Dp < 16 km sec^'kpc"') are more efficient (than the 
corresponding RMs) in portraying the NGC 1300 bar (third to sixth 
rows of Fig.|6ll. 

In Figs. I7I9I we plot the isocontours (dashed lines) corre- 
sponding to the minimum dou with respect to the isocontours of 
NGC 1300 (solid lines), for a sample of the SMs of the Models A, 
B and C, respectively. These samples include characteristic cases 
such as those with Qp = 16 km s"' kpc"' and Sip = 22 km s"' kpc"' 
. Note that the corresponding isocontour curves (of model and 
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Figure 10. The "global" dcH index for all the six contour lines presented in Fig.|6] The left panel corresponds to RMs. while the right one to SMs. These dcH 
values represent globally a model. The plotting styles are as in Fig.|6] 



galaxy) are plotted in the same color. From these figures, it be- 
comes evident the nice behavior of the cIqu index in quantifying 
the comparison between the models and the image of the galaxy. 
In almost all cases low (or high) dcH values imply good (or bad) 
resemblance between the corresponding curves. 

Since the doH values for the various isocontours are by defini- 
tion comparable, we can add all the values of each model in order 
to get a unique value describing globally the model. Such an in- 
dex puts of course equally weights to all isocontours and it could 
be low valued in models with no matching curves. Nevertheless, 
it is reliable in cases with many model curves showing moderate 
resemblance with those of the galactic image. In Fig.[TO]we plot, 
for all RMs and SMs, this "global" dcH index for the isocontours 
ofFig.lH 

The general conclusions we draw by observing Fig.[TO]can be 
summarized as follows: 

• For Sip < 20 km s"' kpc"' values, the RM (left-hand panel of 
Fig.[TOj with the lowest global dgu index are those corresponding 
to Model A. 

• The RMs corresponding to Model B and C improve their be- 
haviour from low Q.p values towards high Sip values (for Model C 
only up to Sip = 24 km s"' kpc"' .) 

• The overall behavior of SMs corresponding to Model A are 
improved (relative to those of the corresponding RMs) only for Q,, 
values around Sip = 22 km s ' kpc ' . 

• The overall behavior of SMs corresponding to Model B are 
improved (relative to those of the corresponding RMs) for Sip < 
23 km s"' kpc"' . 

• The overall behavior of SMs corresponding to Model C are 
improved (relative to those of the corresponding RMs) only for 
Sip < 14 km s"' kpc"' . 



5 DISCUSSION AND CONCLUSIONS 

We studied the stellar response, under the assumption of a single 
pattern speed, for the potential models of the barred-spiral galaxy 
NGC 1300 estimated in Paper I and for a wide range of Sip values. 
The three potential models of Paper I, correspond to three different 
geometries as regards the mass distribution in the third dimension 
(perpendicularly to the equatorial plane). Our goal was to find out. 



which geometries and which values of Sip were able to reproduce 
the fundamental morphological features of NGC 1300, i.e. the bar 
an d the s pirals. We used the method of response models (RM), as 
in IPatsij i2006l) . So, we start with a uniformly populated disk of 
stars moving initially in circular orbits, which are determined in 
the axisynmietric part of the potential. We integrate the orbits of 
the test particles for many pattern rotations under the full potential 
and we get the density maps of the response models as described in 
Sect.m 

The "by eye" comparison between the RM morphologies and 
the image of the galaxy gives useful qualitative information but 
it lacks of objectivity and quantification. For this reason we in- 
troduced a new index which is a generalized modification of the 
Hausdorff distance. This index quantifies the resemblance of two 
isocontour lines and enables the determination of the best matching 
contour lines (of a model) to each one of the six preselected isocon- 
tours of the deprojected K-band image of NGC 1300 (Fig.[T](. This 
new index testifies the general remarks made by eye and further- 
more it describes the quality of the fitting for the individual com- 
ponents of the galaxy. Moreover, we investigated possible improve- 
ments of the response models as regards their resemblance with the 
NGC 1300 morphology, by assigning weights to the contributions 
of the various energy levels. This method follows the general con- 
cept of Schwarzschild's methods, but it is simpler in the sense that 
it requires weights for the energy levels and not for individual or- 
bits. Some of these modified Schwarzschild models (SM) improve 
the similarity of specific features of NGC 1300 in comparison to 
the original response models. 

Below we enumerate our conclusions: 

(1) The spiral structure of NGC 1300 is reproduced: 

• Best by the RM corresponding to the potential of the pure 
2D geometry and for Sip values around 16 km s"' kpc"' . The big 
advantage of this model is that the spirals of RM and galaxy have 
similar pitch angles. The density enhancement of the the RM at 
the spiral region is along the same direction as the corresponding 
arms of NGC 1300. 

• Acceptable by the SM corresponding to the potential of the 
pure 2D geometry and for Sip values around 22 km s"' kpc"' . 
However, the drawback of this case is that the spirals, especially 
the one appearing at the lower part of the figures, have different 
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pitch angles than that of the galaxy. The spiral arms of the model 
share common regions with the NGC 1300 arms but essentially 
cross each other. They also add additional features beyond the 
region within which the NGC 1300 barred-spiral structure ex- 
tends. 

• Moderately by the RM and SM corresponding to the po- 
tential of the 3D cylindrical geometry and for fl,, values around 
22 km s"' kpc"' . However, the remarks about the pitch angles of 
the spirals we mention above for the 2D case, yield for the thick 
disc 3D case as well. On top of that the upper arm is not nicely 
reproduced by this model. 

(2) The ansae-type bar of NGC 1300 is reproduced: 

• Best by the RM and SM corresponding to the potential of 
the pure 2D geometry, for fl,, values around 22 km s"' kpc"' 
(very good is also the da^ index of the RM corresponding to 
the potential of the 3D cylindrical geometry again for fl,, around 

22 km s"' kpc"' . However, the ansae are reproduced in the cor- 
rect regions only in the 2D case) 

• Very well by the SM corresponding to the potential 
of the 3D cylindrical geometry and for Q.,, values around 
17 km s"' kpc"' . 

• Moderately by the SM corresponding to the potential 
of the 3D cylindrical geometry and for Q,, values around 

23 km s"' kpc"' . 

(3) The inner oval-like shape structure of NGC 1300 is repro- 
duced: 

• Best by the RM and SM corresponding to the potential of 
the 3D spheroidal+cylindrical geometry and for almost all £1,, 
values 

• Well by the SM corresponding to the potential of the 3D 
cylindrical geometry and for Q,, values around 22 km s"' kpc"' . 

• Moderately by the SM corresponding to the potential of the 
pure 2D geometry and for f2,, values around 23 km s"' kpc"' . 

All the above indicate the tendencies of the models to improve 
their similarity with the K-band image of the galaxy for a given ge- 
ometry or for a certain £1,, range. However, the relation between 
"good" models and the NGC 1300 morphology is not unambigu- 
ous. It is obvious that some geometries support better certain mor- 
phological features. This may suggest the appropriate geometry for 
the description of specific regions of the galaxy (e.g. the central 
area should be considered as a spheroidal). Another issue is that 
a morphological feature may be reproduced equally well by mod- 
els that diff^er in their pattern speeds. In addition we remind our 
assumption of the time-independency of the potential and a con- 
stant time-independent O,,. The study of the stellar dynamics of 
our system under these assumptions not only helps us see some ba- 
sic trends in the dynamical behaviour that favours the formation of 
the one or the other morphological feature, but after all it will in- 
dicate the necessity of different assumptions. The evolution of the 
models under a time-dependent potential and Q.p is very dubious in 
the response model approach, since it requires specific assumptions 
for the evolution laws for both potential and Sip. The assumptions 
are free and therefore the number of the possible combinations is 
large. Such a study is beyond the scope of this paper and should be 
attacked by means of A'-body simulations. However, in order to get 
a crude impression about how dilferent can be the resulting mor- 
phologies if a basic parameter varies with time, we test a scenario 
with time dependent Q.p in a fixed potential. We have chosen for 
our calculations a pure 2D geometry. The choice is justified by the 



fact that the RM corresponding to the pure 2D geometry portrays 
on the one hand the spiral structure for Q,, » 16 km s"' kpc"' and 
on the other hand the ansae-type bar for fl,, x 22 km s"' kpc"' . 
We adopted a rather arbitrary evolution law for the increase of Q.p, 
which reads 



Slpit) = 19 + 3 tanh -(f - 15) km s^' kpc^' , (6) 



where the time t is measured in periods corresponding to a pat- 
tern rotating with £2,, = 23km s"' kpc"' . The integration time 
of the orbits was up to / = 30. In this numerical experiment the 
initial value of Q.p = 16 km s"' kpc"' , while the final one is 
X 22 km s"' kpc"' . In the left-hand columns of Figs. [6l and [Tol 
we have plotted (with star symbol or light blue lines in the on- 
line version) the doH values corresponding to each snapshot's Q.p 
value. In these figures we observe that for the time corresponding 
to Clp X 17 km s"' kpc"' and fl,, « 22 km s~' kpc~' we get a good 
resemblance with many isocontours of NGC 1300. This could be 
an indication that the present NGC 1300 structure is just a snapshot 
in an evolutionary scenario, according to which the observed mor- 
phology undergoes successive transient phases. However, whatever 
good results we find, we find for increasing Q.p from one value to 
the other. The opposite procedure, i.e. decreasing Q,, according to 
the law in Eq. [6]does not lead to equally nice results. Decrease of 
Qp during an A^- body simulation is frequently observed due to dy- 
namical friction l|Chandrasekhad [ 19431 : bonner & SundeUus|[T993l: 
iDebattista & Sellwooj|l998l) . but there is no obvious physical rea- 
son for the opposite. In our model the particles that participate in 
the spiral structure stay for longer times at the outer regions of the 
disc, i.e in regions that evolve dynamically slower than particles 
at the central region that contribute to the bar structure. So, by in- 
creasing Q.p, the particles that stay longer at the outer disc do not 
have enough time to feel the potential. On the other hand particles 
staying at the central region of the system respond fast during the 
period our simulation reaches the fast Q.p domain. 

Concluding, we remark that taking into account all the models 
we constructed as well as the techniques we used to assess them, 
there are clearly two Clp intervals, where we find morphological 
features of the models resembling the structures of NGC 1300. 
These are Clp values close to 16 and 22 km s"' kpc"' . 

Comparing this result with pattern speeds proposed by 
other authors for NGC 1300 tenglanj[l989l : iLindblad & KristerJ 
1996 ), there is a relative a greement between the fast model 
of ILindblad & KristenI ( Il996h and the group of models around 
n„= 22 km s~' kpc ' , \ yhich we find matching the NGC 1300 
bar (' Lindblad & KristenI j 19961) propose Op=20 km s ' kpc ' for 
a distance of the galaxy D=20 Mpc). The alternative solution pro- 
posed by the same authors (£1,,=12 km s"' kpc"' ) is in a range 
of Sip values, where the bars of our RMs are much thicker than 
the one of NGC 1300 and the stellar spirals are absent. As re- 
gards the work o f lEnglanj ( fT989l) . the proposed pattern speed gives 
equilibrium points in a narrow zone compatible again with the 
Slp=22 km s"' kpc"' values, as it will become evident in Paper 
III. 

In any case, the understanding of the observed NGC 1300 
morphology goes through the understanding of the underlying dy- 
namical mechanisms and this can be done only by studying the 
orbital behavior of the "successful" models. This will allow us also 
the comparison among them, and is done in Paper III in this series. 
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